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^ ' Abstract. We have theoretically studied the effect of dissociative autoionization on the photoelectron energy 

O 

Jh ' spectrum in (1 + 2)-photon above threshold ionization (ATI) of H2 molecules. We have considered excitation 

d ' from the ground X'^Yj'^{v = 0, j) state to the doubly excited autoionizing states of ^I]+ and ^n„ symmetry, 

O ■ via the intermediate resonant -B^ Sj {v = 5,j) states. We have shown that the photoelectron energy spectrum 

^~>' is oscillatory in nature and shows three distinct peaks above the photoelectron energy 0.7 ev. This feature 

■ has been observed in a recent experiment by Rottke et. al., J. Phys. B 30, 4049(1997). 
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I. Introduction 

Studies on dissociative photoionization in molecular hydrogen is being persued for several years, almost two 
decades, both experimentally and theoretically. In this process the molecule is excited above the dissociation 
threshold of the molecular ion by the single or multiphoton absorption and the doubly excited autoionizing 
states play a significant role to determine the photoelectron and the atomic ion spectrum. In most of the 
single photon experiments (Chung et al 1993, Dehmer and Dill 1978, Gardner and Samson 1975, Glass- 
Maujean 1986, Glass-Maujean et al 1995, He et al 1995, Ito et al 1988, 1996, Latimer et al 1992, 1993, 
Strathdee and Browning 1976) the molecule is excited from the ground ro-vibrational level of the ground 
electronic X^S+ state to the dissociation continuum above the ionization threshold leaving the signature of 
autoionization on the dissociative ionization yields. In most of the multiphoton experiments (Anderson et al 
1984, Cornaggia et al 1986, Normand et al 1986, Pratt et al, 1986, 1987, Verschuur et al 1988, 1989, Xu et 
al 1989), the molecule is excited from the ground state to just above the ionization continuum (and below 
the dissociation continuum) or above the dissociation continuum, via intermediate resonant rovibrational 
level of _B^S+, E,F^Y,'^ or C^Hu states. In the above experiments, when the molecule is excited below the 
dissociation threshold (i.e. when only the bound molecular ions are produced in different rovibrational levels), 
the role of doubly excited states on rovibrational branching of the molecular ion has been studied; also, when 
the molecule is excited above the dissociation threshold (i.e. when atomic ions and neutrals are produced 
as dissociative ionization products together with the bound molecular ions), the effect of autoionization 
via doubly excited states on the atomic ion spectrum and the photoelectron spectrum (PES), have been 
studied. Several theoretical studies ( Chupka 1987, Ganguly et al 1986, Ganguly and Rai Dastidar 1988, 
Hazi 1974, Hickman 1987, Kanfer and Shapiro 1983, Kirby et al 1979, 1981, Rai Dastidar et al 1986) have 
revealed that autoionization via doubly excited autoionizing (Al) states play a dominant role in determining 
the autoionization lineshape and the spectrum of ionization yields i.e. photoelectrons and photoions. It has 
been shown that the interference of autoionization with ionization (Fano 1961), leaves its signature on (i) the 
autoionization lineshape, (ii) the rovibrational branching of the molecular ion in the ground state and (iii) 
the spectrum of the dissociative ionization products (i.e. atomic ions and photoelectrons), for the excitation 
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above the dissociation threshold of the molecular ion. 

For the single photon excitation from the ground state, only the ungerade states can be accessed and 
the lowest of these series is situated approximately 25.5 ev above the ground state. In this case the Franck- 
condon transition occurs at very low intenuclear distances (around R = lAa.u.) and at very high continuum 
energy. But in case of resonance enhanced (m+n)-photon transition, the molecule is selectively excited 
to an intermediate rovibrational level of a Rydberg state and hence the Franck-condon trasition to the 
autoionizing state occurs at large internuclear distance. Therefore, in this REMPI ( resonance enhanced 
multiphoton ionization ) technique, the initial condition for dissociative ionization can be changed and the 
molecule can be excited to different continuum energies, much below the energy, reached in the single photon 
transition. Recently, experimental study (Rottke et al 1997) on (3 + 2)-photon dissociative ionization of H2 
molecules has been done by selectively exciting the i3^S+(w = 5, j) levels, from the ground X^S+(w — 0, j) 
states by three photon absorption and then to the dissociative continuum by two-photon absorption from 
these intermediate resonant levels of _B^S+(w = 5, j) state. Angle resolved photoelcctron spectrum (PES) 
and ion spectrum have been studied and it has been found that the PES gives oscillatory structure with 
increase in energy. 

In the present study, an attempt has been made to explain theoretically the origin of such oscillatory be- 
haviour of PES. In this calculation we have considered (1 -I- 2)-photon dissociative ionization of H2 molecules. 
The transition scheme considered here is shown schematically in Fig. 1. By the single photon transition, 
the molecule is first excited resonantly to the i?^S+(w = 5,j) state from the ground X^S+(w = 0,j) state 
of the H2 molecule and by the subsequent two-photon transition the molecule is excited from this inter- 
mediate rovibrational level of -B^S+ state to the dissociative continuum of the molecular ion in the ground 
state (X^S+). In this transition, the doubly excited states of ungerade symmetry are excited and hence the 
dissociative autoionization process via these AI states contributes to the PES. Our aim here is to show that 
the dissociative autoionization process can give rise to oscillatory photoelcctron spectrum, as observed in 
the recent experiment (Rottke et al 1997). 

II. Transition Schematics 



In the experiment, the H2 molecule is first excited from the ground X^S+(f = 0) state to the B^Yi^{v = 5) 
state by three photon absorption considering three transition schemes: 

(z)P(l) : X^^+{v = 0,j = 1) ^ i?iS+(i; = 5,j = 0) 

{ii)R{0) : X^^+{v = 0,j = 0) ^ B^^Uv = 5,j = 1) 

(m)i?(2) : X^^+{v = 0,j = 2) ^ B'S+(i; = 5,j = 3) 

and then by two photon transition to the dissociative continuum and autoionizing states of ungerade 
symmetry. In our calculation we have considered (Fig. 1) single photon transition from the ground X^S+(u = 
0,j ~ 0,1,2) levels to the different rotational levels of i?^E+(u = 5) state by three transition schemes 
mentioned above and then two photon transition from _B^S+(u = 5) state with rotational quantum numbers 
J ~ 0,1 and 3 to the continuum above the dissociation threshold. In this process, excitation to four 
autoionizing states, two lowest AI states of ^E+ symmetry and the two lowest AI states of ^n„ symmetry in 
the Qi series has been considered. To get the photoelectron spectrum, the autoionization decay from these 
AI states to the different dissociation continuua of the ground ^^S+ state of the hydrogen molecular ion 
has been studied. 

III. Theory 

We have used the resolvent operator technique to obtain the dissociative ionization and autoionization rate. 
In this technique a set of equations for the matrix elements of resolvent operator G{z) — —rn, with respect 
to product states, such as \g)\n), \i)\n — 1), \a)\n — 2), |ci)|n — 2), \b)\n — 3} and |c2)|n — 3). Here, \g) 
and \i) are the ground and the intermediate resonant states of the bare molecule, \a) and |6) are the AI 
states of different symmetry excited by the single photon and two photon absorption from the intermediate 
state of the molecule respectively and the continuua |ci) and |c2) are excited by the single photon and two 
photon absorption from the intermediate state respectively (see Fig. 1). Hence the continuua are of different 
symmetry. The continuum |c2) will be of the same symmetry as the intermediate resonant state \i) and 
the continuum |ci) will be of opposite symmetry. The \n) 's are the photon number states. H is the total 
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hamiltonian of the inolcculc+photon system, z is the complex energy of the system and the imaginary part 
of z gives the decay rate of the system. For the (1 + 2)-photon transitions considered here, the equations for 
the matrix elements of resolvent operator are given as follows: 

(Z - Eg)Ggg - Dg^G^g = 1 (1) 

[z — EijGig — DiaGag — I Dic^GcigdEc^ — DigGgg = (2) 

{z - E^)G.g J V..,G.,gdE,, - D^,G.g J2 Da,G,g - J D^^^G^^gdE^^ = (3) 

b 

{z - E,)G,g - D,^G.g - j D,^,G^,gdE^, - / V,.^G.^gdE^^ = (4) 

\Z — Ec-^)Gcig — Dc^iGig — Vc-^aGag — / ^ Dc^bGhg — I DciC2Gc2gdEc2 — (5) 

b '' 

\Z — Ec2)Gc2g ~ Dc2aGag ^ / D c2CiGcigdEc^ — y ^ Vcabtrbg = (6) 

•^ h 

where Eg — Sg + nhuj, Ei = Si + (n — l)hLo, Ea = ^a + {n — 2)huj, _Et, = et, + (n — 3)hLo, Ec-^ ~ e^ + (n — 2)tko 
and Ec2 = Cca + {n — 3)hLo are the energies for the product states, i.e. the ground, intermediate, autoionizing 
states and the two continuua respectively. Here e's are the energies for the bare molecular states. -Djfe's 
are the dipole transition moments between the product states j and k. The summation over b indicates 
that contribution from all the autoionizing states excited by this transition schemes has been considered. 
Formally obtaining Gcig and Gc2g from equations (5) and (6) and substituting in equations (2), (3) and (4) 
one can derive equations as follows: 

^gGgg - DgiGig = 1 (7) 

ZiGig — DigGgg — BiaGag — / ^ CifeGftg = (8) 

b 

ZaGag — BaiGig — 2_^ KabGbg = (9) 

b 

ZbGbg — KbaGag — CbiGig — (10) 

where, 

Zg = z-Eg 

'i 

Z^^ z - Ei - Si + -7i 



1 1 



1 1 

Zi) = z — Eh — Sh-\- -Th — St -\- -75 + 2Rtc2Cib 



^ib — ^ib ^^ ^ic\C2b 
J^ab — ^ab ^acic^h ^ac^cih 



1 f D,,,D^ 



^''-A T^eT"^^^^ 



Qab = Dab + M ^^^^^dE,, +\ I 2li52}^dE, 



and 



AJ z-E,, '-' AJ z-E,, 



_ 7r2 

J^jCpCgk — '~r\j^jCp-L-'CpCg^Cgk\ 



A 

j,k corresponds to bound product states; Cp^q corresponds to product states for continuua; Xc k and Yjc 
correspond to coupling terms between the continuum Cp,q and the bound states designated by k and j 
respectively. Here A = i_l7i-^ I dc c I ^ ' ^i ^^^ ^i ^^^ ^^® photoionization width and ac-Stark shift, divided by 
the factor A, of the bound states designated by j respectively; Ta,b and Sa.b are the autoionization width and 
the shift, divided by the factor A, of the autoionizing states designated by a and b respectively, due to the 
configuration interaction couplings Vaci and Vbc2 with the continuum adjacent to the AI states. The terms 
involving integration over continuum (as in Qia, Qit, Qab etc.) can be expressed as the sum of the principal 
value part and the imaginary part in the usual manner. In the present work we have neglected the principal 
value part and used the pole approximation to derive the final formula for autoionization rate. Equations 
(8) -^ (10) can be rewritten in terms of the detunings {6i, 6a and 6b) from the intermediate and two AI 
states respectively by substituting z = Zg + Eg in the above equations. By solving above equations one 
can formally obtain the matrix elements for the resolvent operator and by inverse Laplace transformation 
of these matrix elements one can get corresponding matrix elements for the evolution operator. Hence the 
probability for ionization can be written as: P{t) = 1 — \Ugg{t)\'^ — \Uig{t)\'^ — \Uag{t)\^ — \Ubg{t)\'^. 



In the weak field limit i.e. considering dipole transition strength from the ground state to the inter- 
mediate state to be much weaker than that for the ionizing transitions and the autoionizing decay, one can 
write down Ggg = „ \^ , where 



(J2 



B..+J: 



b Zfc 
T 



^-E.^ 



and hence the rate for ionization ^ = —2ImZQ. 

By neglecting the higher order contributions, it can be shown that the rate for dissociative autoionization 
via the AI states \b) can be written as: 

where C — \Zi — ^^ ^^|^ ; ^^ can also be expressed as 



dt 



= Mgc,F (12) 



where, F — tt^u-^-ji, is the lineshape factor and the modified dissociative autoionization matrix element 
square A^gc2 = l^'giPlAci l^i'cibPlVbcs |^ 

IV. Calculations 

In the present study we have considered contribution from the dissociative autoionization via two lowest 
AI states of ^Ej symmetry and two lowest AI states of ^n„ symmetry. The energies and autoionization 
widths for these four states are available in the literature Guberman 1983, Kirby et al 1981, Tennyson 
1996). As mentioned before, we have considered here (l+2)-photon transition from the ground state to the 
autoionizing states via the intermediate resonant i?^S+ state. Therefore, to calculate the two-photon term 
via the continuum, one needs to have the electronic wavefunctions for the AI states. The wavefunctions 
for two lowest AI states of ^Ej symmetry have been obtained by minimizing the energy variationally. The 
molecular basis set used is a multi-configuration minimal basis constructed with linear combination of Slater- 
type atomic orbitals. We have used seven configuration wavefunctions constructed with atomic orbitatls of 
s and p symmetry and the principal quantum number upto 4. At each value of R (internuclear separation), 
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all the exponents have been varied to obtain minimum energy. The energies obtained are in good agreement 
with the previous values (Guberman 1983) except at lower values of intcrnuclear seperation R < lAa.u.. 
Calculated energies for the second ^S+ AI state differ by < .006 a.u. from the previous values (Guberman 
1983), except at i? = 1.0 a.u., where the difference is .017 a.u.. For the first ^S+ AI state, the difference in 
energies (our results and Gubcrman's results) are < .0086 a.u., except at _R = 1.0 and i? = 1.4 a.u., where 
the differences are .023 a.u. and .012 a.u. respectively. The electronic wavefunctions for the ground and 
the intermediate states have been obtained from our previous calculations (Khatun and Rai Dastidar 1995). 
For the calculation of the complex nuclear wavefunctions using complex potential for the AI states, we have 
used modified Numerov method (Lane and Gcltman 1967, Rai Dastidar and Rai Dastidar 1979). 

To calculate two-photon dipole transition moment between i?^Sj state and the AI states (i.e. to 
evaluate Qn,; see section III) we have used the pole approximation to integrate over the continuum energy 
(-Bci ) , and hence considered contributions from the continuua associated with the bound vibrational levels of 
the X'^Yi't state of H2 ion, which are energetically below the continuum energy excited by the first photon 
from the -B^S+ state. We have used the multipolar gauge to evaluate dipole transition moments, but one 
can also use the minimal coupling gauge. Since, in the present work, the two-photon term acts as a constant 
(with respect to PE energy) multiplying factor to the autoionization amplitude-square (see equation 12), 
the structure of the PES is independent of the choice of gauge. However, the scaling factor for PES can be 
different, depending on the magnitude of two-photon term in two different gauges. For the present study 
we have calculated the two-photon dipole transition moments for the transition from the -B^S+ state to 
the two lowest AI states of ^S+ symmetry only. It is found that the second AI state of ^S^ symmetry is 
energetically very close to the lowest AI state of ^n„ symmetry (Guberman 1983, Tennyson 1996)) and the 
difference between the values of effective principal quantum numbers for these two states are also very small 
(Tennyson 1996). Therefore for the present calculation we have used the same two-photon dipole transition 
term for the transition to the lowest AI state of ^n„ symmetry as that calculated for the transition to the 
second AI state of ^S+ symmetry. For the second ^n„ state we have used the same two-photon term and 
it is found that the contribution from this second ^n„ state is much smaller than that from the second ^E+ 
state and the lowest ^Hu state. 
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The electronic configuration interaction term between AI state and the continuum has been derived as 
Vbc2{R) — V 2 — ^' ^^"sre autoionization widths Ti,(R) have been obtained from the literature (Tennyson 
1996) and this value of Vf,c2 has been considered to be arising from I = 1 partial wave only. This is a 
valid approximation, since Z = 1 is the most dominant contribution for autoionization from ^S+ AI states 
(Kirby et al 1981). To obtain autoionization amplitude, integration over nuclear coordinates have been 
done considering explicit dependence of the electronic transition amplitudes on the internuclear seperation. 
Similarly, to obtain dipole transition amplitudes, explicit dependence of electronic dipole transition moments 
on the internuclear separation (for bound-bound transition) and also on the photoelectron energy (for bound- 
continuum transition) has been considered for nuclear integration. 

The autoionizing states lying above the dissociation continuum of the molecular ion can decay into the 
dissociation continuua lying between the final continuum energy excited by two-photon transition from the 
i3^S+ state and the dissociation threshold of the ground state (Ar^S+) of the molecular ion. Hence the total 
energy Et, can be expressed as Et ~ E^ + -De, where Eg is the photoelectron energy and Dg is the energy of 
the dissociation continuum of the molecular ion. We have calculated the dissociative autoionization matrix 
element for the transitions to a range of dissociation continuum of the ground state of the H^ ion, chosen 
(arbitrarily) to lie on a grid of .01 ev spacing. 

V. Results and Discussions 

A schematic diagram for the transitions considered here is shown in fig.l. For different molecular potential 
energy curves we refer to Fig.l of Rottke et al, 1997. The two continuua shown in Fig. 1 are of different 
symmetry and correspond to the first ionization threshold (i/j (AT^S^) + e) of the neutral molecule. Results 
for three different transition schemes (-P(l), ^(0) and _R(2)), i.e. for ionizing transitions from J = 0, J = 1 
and J = 3 of i?^S+(u = 5) state are shown in Figs. 2,3 and 4 respectively. In these figures we have plotted 
the dissociative autoionization rate divided by I^F ( where I is the laser intensity and F is the line shape 
factor, see equ. 12), as a function of photoelectron energy. The factor PF, which is independent of the 
photoelectron energy, merely plays the role of a scaling factor as far as the shape and the structure of the 
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PES is concerned. The experimental PES is given in arbitrary units and here we compare our theoretical 
PES (scaled by the factor I^F) with experiment. It is to be mentioned here that, for photoeleetron energies 
less than 0.7 ev, the contribution from the ionization of neutral atoms is dominating (Rottke et al 1997) 
and hence the effect of dissociative autoionization will be prominent only for energies > 0.7 ev. Therefore, 
we will discuss the features obtained in PES for energies above 0.7 ev ( experimental results also show PES 
above 0.7 ev, see fig. 2 of Rottke et al, 1997). 

For this calculation we have taken A ~ 1, since the value of |-DciC2p is much less than unity for this 
calculation. The continuum-contimnim transition amplitude (|£'cic2l) ''^iU be significant only for photoelee- 
tron energies very close to each other (Gordon 1929) and for the present study the photoeleetron energies 
for the transition to the lower continuum ( due to the absorption of single photon from the i?^S+(u = 5) 
state ) are less than .02 a.u. (i.e. less than 0.54 ev ). Hence the value of ji^cicaP will be significant for 
photoeleetron energies less than 0.54 ev for continuum-continuum transition. To get an estimate, we have 
calculated electronic dipole transition moment for continuum-continuum transition (Gordon 1929) and the 
overlap between the three lowest vibrational levels of molecular ion in the ground state, with different dis- 
sociation continuum of the same. It is found that the contribution from this term is orders of magnitude 
smaller than unity (Rai Dastidar 1996). 

In the figures (2,3 and 4), individual contributions to dissociative autoionization from four ( first two 
"'^E+ and first two ^n„ ) autoionizing states have been shown, together with their sum. It is found that the 
PES is oscillatory in nature and the dissociative autoionization from the second Al state of ^S+ symmetry 
is the most important contribution. We have checked that the oscillation is essentially the same as in the 
Franck-Condon overlap (mod-squared) between the complex nuclear wavefunction of the autoionizing state 
and the nuclear continuum of the molecular ion in the ground state, as a function of photoeleetron energy. 
The PES for autoionization from three Al ( second ^S+ and the two ^n„) states show three peaks above the 
photoeleetron energy greater than 0.7 ev and the peak positions are almost the same for these three PES. 
But the contribution from the lowest ^S^ state is much smaller than that from the second ^S+ state and 
shows an out of phase oscillation. The presence of three peaks in PES, above 0.7 ev of photoeleetron energy, 
has also been experimentally observed (Rottke et al 1997). From the figure for PES (Fig. 2 of Rottke et 
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al 1997), the peak positions can be approximately calculated and we find that the peaks are shited towards 
higher energy than that obtained in the present calculation. These shifts are different for three transition 
schemes and the values are approximately within 0.23 — 0.24 ev, 0.17 — 0.2 ev and 0.04 — 0.09 ev for the first, 
second and the third peak respectively. It is to be mentioned here that the shifts of the autoionizing states 
due to autoionization decay have not been considered in our calculation. By using the approximate formula 
(Ganguly and Rai Dastidar 1988, Olsen 1982) Sb — ^^^b, one can calculate the approximate values of 
shift (Sb) for the second Al state of ^S+ symmetry (using the values of width tabulated by Tennyson, 1996). 
Typical values of shift varies from —0.02 to —0.16 ev, in the range of internuclear separation, R = 3.0 to 4.0 
and since this value is subtracted from the energy of the Al state, electronic potential energies will be shifted 
upwards. Hence the peak positions in figs. 2, 3 and 4 will be shifted towards higher energy and this shifting 
will be different for different photoelectron energies. The values of shifts are 0.2, 0.16 and 0.14 ev for the 
first, second and the third peak respectively. It is found that, after shifting, the peak positions are almost 
the same as those observed in the experiment. It is within .03 to .04 ev, .01 to .04 ev and .05 to 0.1 ev for 
the first, second and the third peak respectively. But for comparing our results with the experimental PES, 
the resolution factor for the experiment has to be taken into account (in particular for determining the peak 
positions). It has been mentioned by Rottke et al, 1997 that the resolution is .08 ev at 2 ev photoelectron 
energy but that for energies lower than 2 ev has not been mentioned. We expect that the resolution will be 
better as the photoelectron energy is decreased and we find that the agreement is better for the first two 
peak positions than the third one. 

By comparing figures 2, 3 and 4, it is found that the PES for three transition schemes are of the same 
order of magnitude and the PES for transition from J = 1 level (Fig. 3) is greater in magnitude than that 
from J — (Fig. 2) and J — 3 (Fig. 4) levels of i3^S+(w ~ 5) state. However, in the experiment, although 
the contribution for these three transition schemes are of the same order (in arbitrary unit), photoelectron 
yield for transition from J=0 level of _B^S+(f = 5) state is greater than that from other two transitions. 
By analyzing our results, we find that the (l+2)-photon term for transition to the second ^S+ Al state is 
greater for excitation from j = level of the ground Ar^S+(f = 0) state than that from other j levels. In 
the presnt work, we have done ab-initio calculations to obtain the PES for dissociative autoionization from 
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two ^E+ AI states. But to obtain PES from two ^n„ AI states we have used the same (l+2)-photon term 
as calculated for second ^S+ state. Hence the scaling factor for contribution from these two AI states (^n„) 
may become different for three transition schemes, if the actual values of (l+2)-photon term for excitation 
to these two ^n„ states are obtained. Moreover, the values of (l+2)-photon term may differ in two different 
gauges (minimal coupling and multipolar). Therefore, at this stage, this difference between our results and 
the experiment may arise, as far as the relative magnitudes of three PES (for transitions from J = 0, l,and 
3) are concerned. Nevertheless, further study in this field is warranted. 

In conclusion, we have shown that the oscillatory nature of the experimentally observed PES in (3+2)- 
photon above threshold ionization of H2 molecules, can be accounted for by considering the dissociative 
autoionization of the doubly excited autoionizing states (of ungerade symmetry) above the dissociative 
ionization threshold. 
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Figure Captions 

fig. 1: Schematic diagram for (1 + 2)-photon above threshold ionization. Two contimuia |ci) and |c2) 
are of gerade and ungerade symmetry respectively and correspond to the first ionization threshold 
{H2 [X'^Y,'^) + e) of the hydrogen molecule. 

fig. 2: Dissociative autoionization rate divided by I^F (see text, equ. 12), where I is the laser intensity 
and F is the lineshape factor, has been plotted as a function of photoelectron energy, for two-photo 
ionizing transition from rotational level, J=0 of _B^S+(w = 5) state. Contributions from four AI 

states and their sum have been shown separately. Legends for the curves are: contribution 

from the lowest ^S+ AI state; **** same from second ^S+ AI state; ++++ same from lowest ^n„ 
AI state; <)<0<) same from the second ^n„ state and - sum total of the four contributions. 

fig. 3: Same as in Fig. 2, for transition from J=l level of i3^S+(w = 5) state. 

fig. 4: Same as in Fig. 2, for transition from J=3 level of i3^X!+(u = 5) state. 
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